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Abstract 

In this paper, we use the cosmokinematics approach to study the accelerated expansion of the 
Universe. This is a model independent approach and depends only on the assumption that the 
Universe is homogeneous and isotropic and is described by the ERW metric. We parametrize the 
deceleration parameter, q{z), to constrain the transition redshift (zt) at which the expansion of the 
Universe goes from a decelerating to an accelerating phase. We use three different parametriza- 
tions of q{z) namely, g,(z) = qi+ q 2 Z, g„(z) = q 2 + q^ ln(l -f z) and g,„(z) = ^ + A joint 

analysis of the age of galaxies, strong lensing and supernovae la data indicates that the transition 
redshift is less than unity i.e. zt < 1. We also use a nonparametric approach (EOESS-i-SIMEX) to 
constrain zt- This too gives zt < I which is consistent with the value obtained by the parametric 
approach. 

1 Introduction 

It is now widely accepted that we are living in a Universe which is accelerating at present |[TJ|2l. 
Measurements of distance moduli of Supernovae (SNe la) along with other independent observa¬ 
tions such as results from the Planck satellite (Cosmic Microwave Background (CMB), Baryonic 
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Acoustic Oscillations (BAO)), etc. support this view |l3l|4l|5l. However, the physics behind this 
late time cosmic acceleration is still a mystery. Cosmologists have come up with different models 
to solve this riddle |16||. Further, the discovery of high redshift SNe la provides strong evidence 
that the Universe was decelerating in the past Q. Therefore, it is very important to find the epoch 
of transition from the slowing down to the speeding up of the expansion of the Universe. 

Earlier Shapiro and Turner (2005) used both pai'ametric and non-parametric methods to con¬ 
strain the deceleration parameter, q{z) and transition redshift, zt using SNe la data IH. Melchiorri 
et. al (2007) also considered different dynamical dark energy models to put constraints on zt ||3. 
As expected, these constraints were model dependent. Further in 2007, Rapetti et al. developed a 
new kinematical method to study the dynamical history of the Universe ifTOl . By using the X-ray 
gas mass fraction measurements, they supported the transition of the Universe from a decelerated 
to an accelerated phase. 

Many other authors have used different data sets such as Gamma Ray Bursts (GRBs), Hubble 
Observational Data H{z), BAO, CMB, Galaxy Clusters , lookback time etc. to reconstruct q{z) 
or to constrain the transition redshift |[m[T^[T3l[T4l[T5l[T^[T7l[T8l[T9ll20ll^l22ll2^ . Lima et al. 
(2012) have advocated many techniques that may constrain the transition redshift effectively in 
the future ll24]l . In a detailed study, Del Campo et al proposed three different q{z) parametrizations 
and put constraints on the model parameters by using SNe la, H{z), BAO/CMB data sets ESll . 
Recently, Vargas dos Santos et al. (2015) chose kink-like parametrizations of q{z) and obtained 
constraints on the free model parameters especially the duration of the transition. They also used 
SNe la, BAO, CMB and H(z) data E^ . 

At present there is no well motivated theoretical model of the Universe which can explain 
all the aspects of this accelerated expansion. Hence it is reasonable to use a phenomenological 
approach. In this paper we study a popular cosmographic approach to characterize the properties 
of dark energy through the deceleration parameter. Here we also attempt to reconstruct q{z) 
through a parametric method. This methodology has both advantages and disadvantages. One 
of the plus points of this technique is that it is independent of the matter-energy content of the 
Universe. The only assumption in this approach is that the Universe is homogeneous and isotropic 
on cosmic scales and the FRW metric is sufficient to describe the space-time of the Universe. The 
parametrization technique also helps to optimize future cosmological surveys. On the other hand, 
this formulation does not explain the real cause of the accelerated expansion. Further, the extracted 
value of the present deceleration pai'ameter may be dependent on the assumed form of q{z). 

In this work we use three different parametrizations of q{z) to constrain the transition redshift 
in a model independent manner. We also use information criteria to compare these models. As 
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mentioned above, earlier work based on q{z) parametrizations used SNe, H{z), BAO, CMB etc. 
to study either the model parameters or the transition redshift. In this paper we highlight the 
importance of strong lensing data to study the “dynamic phase transition” in the universe when 
it switched from a decelerating to an accelerating phase. We further supplement this dataset with 
the age of galaxies and the recent data of SNe la (Joint Light Curve Analysis). 

To check the consistency of our results, we also constrain zt using a nonparametric approach. 
Many nonparametric approaches like Principal Component Analysis (PCA), Gaussian Process 
(GP), Non Parametric Smoothing method (NPS) etc have been used in the literature to constrain 
cosmological parameters ETll . Here, we use the LOcally wEighted Scatterplot Smoothing method 
(LOESS). This method locates a smooth curve among the data points without having any informa¬ 
tion about the priors i.e. the functional form estimation is done using raw data only. Eor details of 
the method see Il28l 129113011 . This method does not consider the contribution of the observational 
errors. Therefore, in order to include the effects of errors in the analysis, LOESS is combined 
with SIMEX (SIMulation Extrapolation method) OOll . 

The outline of the paper is as follows. In Section 2, we describe the parametric approach to 
put constraints on deceleration parameter {q{z)) and transition redshift (zt). In Section 3 we adopt 
a nonparametric approach to constrain zt. Discussion is presented in Section 4. Einally, Section 5 
contains some concluding remai'ks. 


2 Parametric Method 


2.1 The Parametrizations oiq{z) 


The reconstructed form of the equation of state of dark energy (codeiz)) is usually written as OTll 



( 1 ) 


n=0 


By analogy, we write the parametrized form of the deceleration parameter as 



(2) 


n=0 


Eor different Xn{z) we get different parametrizations of q{z): 


Redshift: Xn{z) = z” 

Logarithmic: Xn{z) = [ln(l + z)]^ 
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Scale factor: 



2n 


We use three different parametrizations of q{z) corresponding to the above Xn{z) ifTTlI^I^ . 

q^{z)=qi+q2Z (3) 


quiz) = qi + Q' 4 ln(l + z) 


(4) 


= 2 + (TT^ 


(5) 


The first parametrization is the Taylor series expansion of deceleration parameter around z = 0 
with qi as the present value and q 2 as the first derivative of q{z) w.r.t. z. The second parametriza¬ 
tion gives the present value of the deceleration parameter as <73 at z = 0 while at higher redshifts 
it diverges. The third parametrization converges to 4 as 2 : becomes large. 

The Hubble parameter, H(z) can be written in terms of q{z) as 


H(z) = Hq exp 


2 : 



1 + q{x) 
1 + X 


dx 


( 6 ) 


Lo J 

where Hq is the present value of the Hubble constant. 

The corresponding Hubble parameter for the above mentioned parametrizations can be ex¬ 
pressed as- 


Hi{z,p) = Hoil + 2 ;)^+''* ® exp{q 2 z) 


Hniz,p) = Ho{l -I- 2 )^+® exp ^[ln(l -f 2 ;)] 


12 


Hiu{z,p) = Ho{l + z )2 exp 


^5(2:^ + 22;) 


2(1+ 2)2 

Here p represents the model parameters, i.e. the various q'-s. 


(V) 

( 8 ) 

( 9 ) 











2.2 Datasets and Method 


2.2.1 Age of galaxies 

We use two datasets of age of galaxies: Sample A and Sample B. Sample A consists of the age 
of 32 passively evolving galaxies in the redshift range 0.117 < z < 1.845. This dataset is the 
collection of three subsamples |[34ll . 

The first subsample contains 10 early type field galaxies. The SPEED model is used to derive 
the age of these galaxies. It is assumed that stars in an elliptical galaxy have a single metallicity 
value which is an over-simplification. Eurther, the formation of stars in the single burst model 
may bias the age result as the formation of stars takes place at low redshift also. This problem 
can be rectified by considering only the old galaxies in each redshift bin. Jimenez et al. (2004) 
concluded from their study that the systematic errors associated with age of galaxies is not more 
than 10-15% even with these assumptions |[35l . 

The second subsample consists of 20 old passive galaxies released by GDDS (Gemini Deep 
Deep Survey). The spectra obtained in this survey are of very high quality. Eor this subsample, 
McCarthy et al computed the synthetic spectra using the PEGASE.2 model and then compared 
the obtained spectra with the observed spectral energy distribution 1361. Simon et al. further 
reanalyzed this subsample using the SPEED model and obtained the age within 0.1 Gyr of the 
estimate of the GDDS collaboration. Einally, they emphasized that systematics are not a matter of 
concern for this subsample 13^ . Dust depletion may also be a cause of systematic error, but this 
effect is again negligible. 

The third subsample of sample A consists of two radio galaxies 53W091 and 53W069 13711381 . 
While deriving the age for these galaxies it was assumed that the large elliptical galaxy was formed 
in a single burst after which there was no star formation in it. This can cause over-estimation of 
age. However, the effect of the metallicity on the calculated age was found to be negligible. 
Eollowing Samushia et al., we also use 12% error in all the observed age data points l39l . 

We also work with another recent low redshift age dataset (Sample B) which consists of ERGs 
(Euminous Red Galaxies) selected from the SDSS DRV. This dataset consists of 12 points in the 
redshift range 0.05 < z < 0.4 iOl. The full spectrum fitting method is used to obtain the 
age of these galaxies by using the single population synthesis model (GalexEV/SteEib). The 
reason for using the combined spectrum is to improve the signal to noise ratio and to reduce 
contamination. Eiu et al argue that the recent star formation in the galactic center of massive 
galaxies may systematically underestimate the obtained mean age at z < 0.4 HOl . But this bias 
is not very significant for galaxies with a large velocity dispersion. Eollowing this criterion, we 
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also prefer to use subsample IV of the data (Liu G., private communication, 2015). 

We use the minimum chi-square technique to find the best fit model parameters. 

where ni is the number of data points used in the sample, t^^{zi,p) is the theoretical age of 
the galaxy, r is the delay factor or incubation time, is the measured current age of the galaxy 
and (Tj is the uncertainty in the observed age. 

The theoretical expressions for the ages of galaxies corresponding to the three 

parametrizations are as follows: 


1 

l + Z 


ti{z,p) = J 


® exp 


-Qi 


1 — X 


dx 


( 10 ) 


1 

1 + 2 


tuiz,p) = j 


X® exp 


^4 

In 


2 ‘ 

2 


VxyJ 



dx 


( 11 ) 


1 

l + Z 


tui{z,p) = iTf 


-'Jj 


exp 




dx 


( 12 ) 


0 

Here p represents the model parameters, i.e. the various q's. 

In an earlier work, marginalization was done over the delay factor only and Ho was kept fixed 
for fhe kinematic approach ifTTl . The recenf value of Hq determined by SNe and Mega-masers 
and also fhrough a combination of dafa including fhe age of old galaxies af intermediate redshift 
are 73.8 ± 2.4 Kms“^Mpc“^ and 74.1 ± 2.2 Kms~^Mpc~^ respectively Il4ni42]| . However, the 
Planck results suggest Hq = 67.8 ± 0.9 B^ . Therefore instead of choosing fixed 

value of Hq, we minimized w.r.f. Hq in order fo fake care care of fhis discrepancy. We find 
fhaf minimization and marginalizafion are equivalenf in our case. 


= A + Ct^ + 2Bt 


(13) 


So we firsf minimize chi-square over fhe delay facfor(r). The minimum value of chi-square is 

X^ = A-^ (14) 

where 
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is again minimized over Hq 


dA 2B dB 


X = 
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and 


dx^ _ _ 

dHo dHo C dHo 

dA -2D 


(15) 
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dB 
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Putting the values of 


where 


and in Equation (fTSl) and equating to 0, we get 
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The final chi-square obtained after minimizing over both the delay factor and the Hubble 
parameter is given by 


C C{E^ - CD) 


2 . 2.2 Strong tensing 


When light passes near matter, it bends due to the gravitational field. The same phenomenon 
occurs if there is a galaxy or galaxy cluster (lens) in the path of light coming from a bright object 
(quasar) usually termed as the source. If the source, lens and the observer are placed in such a 
way that all three lie in the same line of sight then a ring-like structure, called an Einstein ring, is 
formed and the phenomenon is called strong lensing (SE). Two or more images can form in strong 
lensing (for details see ll44]l l. 

The phenomenon of strong lensing can be used to constrain cosmological parameters. If we 
assume that the Singular Isothermal Sphere (SIS) or Singular Isothermal Ellipsoid (SIE) model 
can be used to represent the gravitational lens, then the radius of the Einstein ring is | 


6 e = dvr 


^sis d^ls 

c2 D„, 


(17) 


Here is the radius of the Einstein ring. Dig is the angular diameter distance between the 
lens and the source, Dos is the angular diameter distance between the observer and the source and 
(7sis is the velocity dispersion of the mass distribution of the lens. 
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We define the ratio of the measured angular diameter distances, i.e. Dig and Dos, as 


'jyobs _ 


D. 


c^9e 


( 18 ) 


Dos ^TTCJ^j, 

We select a subsample from the Sloan Lens ACS (SLACS) and Lens Structure and Dynamics 
survey (LSD). This subsample consists of 70 data points in which lensing occurs due to galaxies 
while for the rest 10 lensing occurs due to clusters. In treating lensing due to clusters, several 
assumptions are commonly made like the whole gas inside is at same temperature, the gas pressure 
and gravity of the relaxed cluster balance each other according to hydrostatic equilibrium and the 
cluster is spherically symmetric. These assumptions may over-simplify the results and hence we 
do not include cluster data in our analysis. We also do not consider the four-image lens systems 
as the SIS lens model produces two images only. Finally we are left with a subsample which 
consists of 36 data points. It is important to note that though the range 0 < < 1 is the 

only physically meaningful one as Dig should always be smaller than Dos, we have also included 
points with T>°^^ > 1 since this is within Icr. The redshift ranges of the lens and the source are 
0.106 < zi < 1.004 and 0.1965 < Zg < 3.9 respectively (for details see ll43]l l. 

To constrain the best fit model parameter, chi-square is written as 


where n 2 = 36 and a-jjobs is the uncertainty in 

The corresponding theoretical quantity can be written as 

■\th 


(19) 


V^^{zi,Zs]p) = 


Dj^ _ Da{zi,Zs]p) 


Da{0,Zs;p) 

In flat Friedmann-Robertson-Walker (FRW) cosmology, the angular diameter distance is 

1 c f dx 


( 20 ) 


Da{zi,Zs-,p) = 


Da{0,Zs;p) = 


{l + Zs)HoJ E{x;p) 


( 21 ) 


2 s 

/ 


dx 


(1 -h Zs) Hq J E{x]p) 
0 


E{x]p) = 


H{x\p) 


Hn 


( 22 ) 


(23) 


In this method, the cosmological model enters thr'ough the ratio of the two angular diameter 
distances and so it is completely independent of the Hubble constant. 
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2.2.3 Supernova la 


We use the latest Supernova la data which consists of 740 data points B9l . This includes data 
from low redshift (z < 0.1), SDSS-II (0.05 < 2 : < 0.4) and SNLS 3 year (0.2 < 2 < 1.0) 
supernova samples. It is assumed that the supernova will have the same intrinsic luminosity at 
all redshifts provided they have identical colours, the same galactic environment and shape. This 
assumption is expressed through the distance estimator which can be written as |[50l 


fJ-obsici, /3; M) = m*^ — M + aXi—(3C 


Here m*^ is the rest-frame B band peak magnitude, and a, (5 and M are the nuisance parameters. 
Xi and C describe the stretching of the light curve and colour at maximum brightness respectively. 
The chi-square for the supernova analysis is written as 



where 713 = 740, pth represents the model parameters and pobs are the nuisance pai'ameters 
present in Pobs^ i-S- « and [5. We fix the values of a and /3 according to the ACDM model as 0.14 
and 3.14 respectively while we minimize over the parameters M and Hq ll50l . 

Chi-square obtained after minimization over M and Hq is 


Xj LA^PthiPobs) — ^ 





Pth = 5(0510 (di) + //0 


where di is the Hubble-free luminosity distance and 


/To = 5(0510 (c/i7o) + 25 


The combined chi-square is the sum of chi-squares obtained from all three datasets i.e. 
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(24) 


2 _ -2 
^tot X.age 


+ XSL + XjLA 


2.3 Results 

We parameterize the deceleration parameter with a two-parameter function iqi{z), qii{z)) and a 
one-parameter function {qni{z)) to understand the expansion history of the Universe.. By fitting 
these functions to the combined set of Age, SL and SNe la data, we get the constraints described 
below. 

• Parametrization I qi{z) = qi + q2Z 

Figure [T] shows 68.3%, 95.4% and 99.7% confidence level ellipses in fhe model paramefer 
plane for fhe combined (Age-i-SL-i-JLA) dafasef. Figure|2]shows fhe evolution of q{z) w.r.f. z wifh 
Ifj, 2fT and 3a errors. 

The besf til parameter values for fhis paramefrizafion are fabulafed in Table [T] The value of zt 
from fhe sfrong leasing dafa is relatively large buf combination of fhe fhree dafasefs (Age-i-SL-i-JLA) 
lowers if fo a value which mafches wifh ofher observations. The value of qi{0) for fhe combined 
analysis is —0.52. 


Table 1: Values of best fit parameters for ( 7 ,( 2 ;) 


Dateset 

Age (Sample A) 

Age (Sample B) 

SL 

JLA 

Age(A-i-B)-i-SL-i-JLA 

Qi 

- 1.68 

-0.052 

-0.85 

-0.57 

-0.52 

Qi 

1.55 

0.54 

0.23 

0.74 

0.53 

9,(0) 

- 1.68 

-0.052 

-0.85 

-0.57 

-0.52 


12.28 

8.60 

64.71 

624.82 

842.31 

xl 

0.41 

0.86 

1.9 

0.85 

1.03 

Zt 

1.08 

0.10 

3.70 

0.77 

0.98 

FoM 

0.71 

0.007 

4.62 

84.52 

99.75 
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Figure 1: la, 2a and 3a confidence level contours of chi-square in the qi — q 2 plane for the combined datasets 
(Age-i-SL-i-JLA). The black dot shows the best fit value of the parameters. 



Figure 2: The variation of qi{z) vs 2 ; with la, 2a and 3a confidence regions. The thick black line is the best 
fit line. The horizontal line represents q{z) = 0. 

• Parametrization II: ^ 11 ( 2 ;) = qi + q^ ln(l + z) 

For the combined dataset, the la, 2a and 3a confidence contours in the q^ — q^ plane are 
plotted in Figure [3l The value of zt obtained from the strong lensing data is once again on the 
high side, but by combining the strong lensing, age and JLA data, we get a value of zt which is 
consistent with other observations such as H{z) etc. within 3a level. The value of (7ii(0) for the 
combined dataset is —0.56. The best fit values of the model parameters are given in Table|2] 
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Table 2: Values of best fit parameters for g„(z) 


Dataset 

Age (Sample A) 

Age (Sample B) 

SL 

JLA 

Age(A-i-B)-i-SL-i-JLA 

Qi 

-1.88 

-0.04 

-0.84 

-0.61 

-0.56 

Q 4 

2.54 

0.54 

0.30 

1.02 

0.83 

<Jn(0) 

-1.88 

-0.04 

-0.84 

-0.61 

-0.56 


12.35 

8.60 

64.76 

624.49 

841.38 

X 

to 

0.41 

0.86 

1.9 

0.85 

1.03 


1.09 

0.08 

15.44 

0.82 

0.96 

FoM 

0.39 

0.006 

2.58 

65.32 

73.90 


15 

10 

L5 

C7 

to 

15 

10 

>13 >1.0 -IS OJ 

q3 

Figure 3: la, 2a and Su confidence level ellipse for the combined observations. Black dot represents the best 
fit value of the model parameters. 

• Parametrization III: ^111(2) = ^ + (T+^ 

The value of best fit parameters for different datasets are tabulated in Table |3] g'iii(O) for the 
joint analysis is —0.8. Figures [5]and[2show the variation of the chi-square with model parameters 
and the evolution of the deceleration parameter with redshift respectively. 
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Figure 4: Evolution of qn{z) with z. The filled regions are the \a, 2a and 3 it confidence level regions 
respectively. The thick black line is the best fit line. The horizontal line represents q[z) = 0. 


Table 3: Values of best fit parameters for qm^z) 


Dataset 

Age (Sample A) 

Age (Sample B) 

SL 

JLA 

Age(A-i-B)-i-SL-i-JLA 

Qs 

-2.58 

-0.57 

-2.14 

-1.3 

-1.30 

QmiO) 

-2.08 

-0.07 

-1.64 

-0.8 

-0.80 


12.47 

8.67 

66.13 

625.23 

841.62 

xl 

0.42 

0.86 

1.9 

0.85 

1.03 

Zt 

1.27 

0.07 

1.07 

0.61 

0.60 


0.5 



0.1 1.5 1.0 LS 10 IS 


Z 

Figure 6: The evolution of qjji{z) with The la, 2a and Su confidence regions are shown by fhe filled 
regions. The fhick black line is fhe besf fif line. The horizonfal line is fhe line corresponding fo q{z) = 0 
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Figure 5: Chi-square vs model parameter variation is shown by the Dot-Dashed line for the combined obser¬ 
vations. Black horizontal lines are the Icr, 2a and 3a lines. 

2.4 Information Criteria 

Information criteria are used to rank different models that are available to describe a given dataset. 

We use the Akaike Information Criterion (AIC) and the Bayesian Information Criterion (BIC). For 
details see 12111521153 



AIC = x^ + 2P 


BIC = + 2Pln(i 

where P and d are the number of model parameters and data points in the dataset respectively. 
X^ is the minimum value of chi square. If d/P < 40, it is good to use the corrected Akaike 
Information Criterion (AICc) 123 ■ 


AICc = AIC -h 


2P{P + 1) 


(d-P-l) 

In our work d/P > 40 and so we calculate AIC for all three parametrizations. We also calculate 
BIC. To compare the models, we calculate the differences given by AAIC and ABIC. 


AAIC = AICj - AlCm 
ABIC = BICj - BICm 

AlCm and BICm are the minimum value of AIC and BIC respectively. If more than one model 
fits the observations equally, then the value of AIC and BIC will be minimum for the model hav¬ 
ing the minimum number of model parameters. 
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• Summary of Information Criteria Results 


The results for different models using the information criteria is given in Table ID 

Table 4: Information Criteria results for Age, SL and JLA data 


Model 

P 

ABIC 

AAIC 


2 

14.11 

4.70 

Quiz) 

2 

13.18 

3.77 

Qiiiiz) 

1 

0.0 

0.0 


• We also use the Figure of Merit (FoM) tool to quantify the robustness of the constraints 
obtained. For two parameters, FoM is defined as the inverse of 95% confidence limit area. Tighter 
constraints in the parametric space give large values of FoM’s. 


3 Nonpar ametric Method 

3.1 Hubble Data 

For this, we use recent Hubble data which includes 30 data points. The redshifts of the passively 
evolving galaxies are known with high accuracy. So these galaxies can be used as chr'onometers 
and can provide the direct measurement of the H{z). This is known as the differential age ap¬ 
proach and 23 data points in this dataset are measured using this approach. The second way to 
measure H(z) is by the clustering of the galaxies. In this method, H(z) measurement is done by 
using the peak position of the BAO in the radial direction as the standard ruler. This approach is 
called clustering. Seven measurements of H{z) are done by clustering. The complete dataset is 
taken from Yun Chen et al. 14^ 


3.2 Method 

Non parametric method includes the study of the neighbourhood points of the focal point. For 
this we need to calculate the smoothing parameter or span (s). This parameter tells us about 
the proportion of observations to be used in the local regression, s ranges from 0 to 1. To find 
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the appropriate value of s, we use the cross-validation method 1301. In this method, the cross- 
validation function {CV) is defined as 

CV{s) = - y^{h{zi) - h{z.i)f (25) 

i 

where n is the number of data points used in LOESS, h{zi) are the data points of the original data, 
h{z-i) are the fitted values obtained after omitting the observation from the local regression 
method at the point of interest i.e. Zi^. CV is calculated for different s i.e. from 0.2 to 1. The 
value of s which minimizes CV is the best smoothing parameter for the particular data set. For 
the data set we use, s = 0.92 is the best value. 

The number of data points that we will use in LOESS can be calculated using the relation 

n = N.s (26) 

where N is the total number of data points. In the data set we use, N = 30 and s = 0.92 so 
n = 27. 

The next step of the method is the estimation of kernel in such a way that the points closer to 
the focal point (zi^) get more weight then the farther points. This is because the points which are 
closer to the focal point are more correlated and so will affect the estimation of h{zi) more than 
the farther points. The weight or kernel function can be calculated using the following tricube 
function 


where 


Wij = 


0 


for \uij\ < 1 
otherwise 


^ 


(27) 


(28) 


A is the maximum difference between the focal point and the last (j*^) element of its window 

i.e. max|zj — Zifl\. 


To calculate the value of reconstructed h{z) i.e. h{zi), we define chi-square as 

n 

= '^ Wij{h{zi) -a-{b* Zi)f (29) 

i=l 

The value of [o -|- (6 * Zi)] which minimizes chi-square gives the value of h{zi). This process 
is repeated for all data points to find the reconstructed value corresponding to all redshifts given 
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in the real data. 


3.3 Results 

The CV versus s plot is shown in Figure |7] This shows that the smoothing parameter value which 
minimizes the cross validation function is 0.92. 
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Figure 7: C V vs s plot for the Hubble data 

Figure shows the variation of H{z)/{1 + z) with z. The scattered points are the values of 
ff(z)/(l + z) from the data with error bars.The blue line represents the reconstructed values of 
H{z)/{1 + z) using the LOESS method. In the LOESS method observational errors are not taken 
into account. This can be done using the SIMEX. The results of using LOESS with SIMEX is 
shown by the orange line in Eigure[8j The detail of SIMEX are given in [SOl. It is clear from the 
plot that the transition from decelerated expansion to accelerated expansion occurs ni zt ~ 0.7. 
Eor reference we plot the H{z)/{1 + z) for ACDM model (green line) which gives zt ~ 0.6. 
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4 Discussion 

In the present analysis we use a sample of age data (passively evolving galaxies+LRGs) which has 
44 points and is therefore larger than the data used earlier. As pointed out earlier, Jimenez et al. 
(2004) showed that the uncertainties associated with the age data may not be more than 10 — 15% 
l|35l. Recently, Jun Jie Wei et al. analysed the data with 24% uncertainty in the observed age of 
passively evolving galaxies in order to take care of some of the unknown systematics which are 
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Figure 8: Vaiiation of H{z)/{1 + z) with z\ the blue line corresponds to the reconstructed values through 
LOESS, the orange line shows reconstructed values from LOESS+SIMEX and the green line is for ACDM 
model. 


not sufficiently understood |[55l . As expected, the constraints obtained on the model parameters 
were weaker. Wei et al. also relaxed the assumption of a uniform delay factor for all galaxies. 
Instead they considered the case in which the delay factor is distributed through a Gaussian and 
a Top Hat distribution. They conclude that the delay factor distribution does not affect the result 
significantly. In the present work, we minimize both the ejfect of delay factor and Hubble param¬ 
eter analytically, something we believe has not been addressed completely in the past. 

Strong gravitational lensing is another important tool used in this work to constrain the tran¬ 
sition redshift. Recent work done by different groups has made this probe more popular to con¬ 
strain the cosmological parameters BSl l56l l57l l58]l . But there are many systematic issues with 
this method. The three important uncertainties associated with this method are the unknown mass 
distribution of the galaxies, stellar velocity dispersion and the line of sight mass contamination. 
Eurther, the data for lensing is drawn from different surveys which may also cause a systematic 
error. 

As the true mass profile of the galaxies is not known, various models have been proposed in 
the literature. But due to the simplicity, most frequently used are the SIS (Singular Isothermal 
Sphere) and SIE (Singular Isothermal Ellipsoid) models. The advantage of considering the SIS 
model is that, to a first approximation, it describes the average properties of the galaxies quiet 
accurately. However, real galaxies show some deviation from the assumed symmetry. Models in¬ 
cluding ellipticity provide a better description of mass distribution but they need a comparatively 
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large number of parameters. 


As mentioned earlier, the total density distribution of the lenses is described by the SIS model 
(stellar and dark matter halo). In principle, the velocity dispersion of the stellar component might 
be different from that of the dark component. This fact is supported by X-ray properties of late 
elliptical galaxies. Treu et al. (2006), after analysing the large sample of lenses of SLACS survey 
ll59l . noticed that the average ratio of cjo and asis (for the total mass) is very close to unity. This 
result suggests that some unknown nexus exists between the stellar and dark matter halo which is 
sometime referred to as the ’’bulge halo conspiracy”. In practice, cJo and asis are related as 

<^SIS = /eCTO 

Here Je is a free parameter which incorporates the uncertainties like the error introduced 
in 6*^; due to the SIS assumption, the rms error between the observed velocity dispersion and 
model velocity dispersion etc. In literature, /e has been treated in a variety of ways. Melia et 
al. for instance, kept it as a free parameter |[57l . Wang and Xu present a new method which 
completely eliminates the uncertainty appearing due to Je II 6 OII . Van de Ven et al have found that 
< Ie < \/T2 |[6n . In this work we assume fE = ^- 

The third systematic is due to line of sight mass contamination |[5^ l62l l6^ . In the present 
work, lenses are assumed to be isolated. In principle the line of sight density fluctuations can 
affect the lens model also. Bar-kana showed that large scale structure (LSS) changes the angular 
diameter distance at a given redshift ll64]| . This shows that if the effect of LSS, is not included 
while reconstructing the lens, it may lead to incorrect conclusions. This effect dominates when 
the source and the lens are at large redshifts. Keeton et al claim that external shear perturbation is 
small due to LSS so its effect on gravitational lenses can be neglected Il65l . A larger dataset for 
strong lensing may help to get rid of this effect. 

The total uncertainty present in V°^^{aQ, Je^Ge) calculated by using the propagation equa¬ 
tion. The uncertainty both in the Einstein radius 9e and the velocity dispersion cJo are ~ 5% 
each |[66l . The error in Je is nearly 6% ll57l 15^16X1 . So the total error in comes out to be 
approximately 16% . 
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5 Conclusions 


More than four decades ago, Allan Sandage (1970) argued that we need a precise measurement 
of the Hubble constant (Hq) and the deceleration parameter q{z) to distinguish between various 
cosmological models 1671. But with the progress in the cosmological observations, it is believed 
that we may require more than two numbers to understand the expansion history of the universe 
l68l . In this regard, the reconstruction of q{z) and the study of the transition redshift are very 
important in describing the evolution history of the universe. 

One of the objective of this work is to include strong lensing (SL) data for the reconstruction 
of q{z) through a parametric approach. There are many advantages of using SL data in the form 
of distance ratio. Firstly, this ratio is independent of the Hubble parameter. Secondly, this data 
can lift the degeneracies in the determination of the constraints obtained from the other observa¬ 
tions. Finally, this method does not depend upon the evolutionary effect of the sources and is also 
immune to dust absorption. In contrast to the SL data, the constraints obtained from the age data 
depend upon the Hubble parameter. Therefore we have minimized the w.r.t. Hq and the delay 
factor analytically, which was not addressed completely in the past ITU . 

The constraints obtained on the transition redshift by using strong lensing, age of the galaxies 
and recent SNe la data through a parametric approach are as follows: 

1. All the three q{z) parametrizations are consistent with an accelerating universe at the 
present epoch and the best fit value of zt is less than 1. 

2. Within 3a confidence level, bofh fhe presenf values of fhe decelerafion parameter and fhe 
fransifion redshifl are in concordance wifh fhe consfrainfs obfained from fhe ofher observafions 
such as H(z), BAO, Clusters efc. 

3. The evolufion of fhe deceleration paramefer w.r.f. fhe redshifl shows lhal Ihere is no slowing 
down of fhe presenf cosmic acceleration as reporfed earlier llti^lTOl . 

4. We observe fhe tension befween fhe resulls obfained from fhe high redshifl ( Sample A) and 
fhe low redshifl age dala ( Sample B). This inconsislency disappears when if is combined wifh 
ofher dafasefs. 

5. The presenf value of q(z) and zt obfained from Ihese paramelrizalions (qi, qjj) have large 
error bars al high redshifl. This may be due lo lack of dala poinfs and more imporlanlly fhe uncer- 
lainfies associafed wifh fhe lensing dala sef. The presenf slrong lensing dala is nol precise enough 
lo conslrain fhe model paramefers lighlly. Therefore, we use fhe combination of age, slrong lens- 
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ing and SNe la data. These observations have different restrictive powers in the parametric space. 
Hence, it is always better to use complementary tools to put constraints on the model parameter. 
This also helps to break the degeneracy inherent in the parametric space. With the addition of age 
data and SNe la data, the constraints on the parameters improve and become consistent with the 
ACDM model. 

6 . We use two information criteria (AIC & BIC) for selecting the best parameterization. ABIC 
greater than 6 is considered as strong evidence against the model with the higher BIC |[5^ . 
The rule remains same for AAIC. In our calculation, AAIC and ABIC are minimum for the 
third parametrization as expected. The value of ABIC is greater than 6 for the first and second 
parametrizations, which indicates a positive evidence against these models. We believe that the 
IC analysis favours the third parametrization. 

In order to obtain the effective constraints using SL data, one must control the associated 
systematics. One very important parameter linked with SL measurements is the observed velocity 
dispersion (Dig/Dos a^). As reported by Schwab et al. 11711 . the actual velocity dispersion is 
luminosity weighted along the line of sight and also over the spectrometer aperture. In addition 
to all these factors, one should also take care of the radial and tangential component of velocity 
dispersion. After including the above factors in the analysis, is replaced by f/^ 


dvrcr, 


0 

This modified rafio fakes care of fhe anisotropy present in the velocity dispersion as well as 
other lens parameters. 

As 2?°^^ is dependent on cjo, A^(a,/3, 7 , 5) and 9e, it is crucial to check the sensitivity of 
these parameters on the If we substitute the mean values of these parameters, we find fhat 
the value of 2?°^^ gets reduced by around 12%. What is more surprising is that this new data 
set of reduced leads to a universe which is decelerating at present, something which is not 
supported by observations. This leads us to conclude that our understanding of the lens quantities 
needs to improve further. 

The point we want to highlight is that the determination of lens parameters as well as the 
observed quantities in the SL data are very crucial. Small changes in the lens parameters may 
change the result completely. Therefore a better understanding of the assumptions and the sys¬ 
tematics associated with the lens pai'ameters must be addressed. 

The obtained constraints are limited by the quantity and the quality of the data for strong lens- 
ing. Nevertheless, the obtained results are meant to highlight the importance of strong lensing as 
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well as the updated age data which gives an independent estimate of the cosmological parame¬ 
ters. At present the known gravitational lenses are of the order ~ 10^. It is expected that next 
generation space telescopes such as Euclid will observe ~ 10^ strong galaxy-galaxy gravitational 
lenses. So in the near future, we can expect that the enhanced data sets will allow us to put tighter 
constraints on the cosmological parameters f/^ . Similarly it is also expected that future surveys 
such as BOSS will observe a large number of quiescent massive LRGs and hence more accurate 
age-redshift data of galaxies will be obtained |[74l . This may further reduce the systematic uncer¬ 
tainties associated with age data. 

The constraints obtained from the nonparametric approach are summarized as follows: 

7. We have applied a nonparametric method (LOESS) to constrain zt by using the Hubble data 
only because this technique can reconstruct the global trend of the observed quantities without 
assuming any prior or cosmological model. So by studying the variation of H{z)/{1 + z) w.r.t. 
z through LOESS method one can constrain zt easily. This method gives a smooth curve using 
a nonparametric regression thr'ough the real scattered data points. The value of zt comes out 
to be about 0.7. This is in agreement with the value we obtained using the model independent 
(cosmokinematics) approach. Eurther, the transition redshift obtained from the LOESS-i-SIMEX 
analysis does not change significantly but it scales down the variation of H(z) /(1+z) as compared 
to its values obtained from the LOESS only. But we did not use Hubble data in the parametric 
approach as this dataset has been used recently for constraining the deceleration parameter 1261. 

The most important point is that knowledge of the nature of the dark sector of the Universe is 
still woefully inadequate. We believe therefore that we should use model independent techniques 
(both parametric as well as nonparametric) with all available datasets for extracting as much in¬ 
formation as possible about the Universe and its structure. 

Acknowledgments 

Authors are thankful to the anonymous referee for useful comments. NR acknowledges fi¬ 
nancial support from the UGC Non-NET scheme (Govt, of India) and the facilities provided at 
lUCAA Resource Centre, Delhi University. NR also thanks Remya Nair (lUCAA, Pune) and San- 
jeev Kumar (DU) for helpful discussions. One of the author (DJ) is thankful to CTP (Jamia Milia 
Islamia, New Delhi) for research support. AM thanks Research Council, University of Delhi, 
Delhi for providing support under R & D scheme 2014-15. Authors are greatful to Lixin Xu, 
Eulvio Melia, Marek Biesiada, Ariadna Montiel and Vincenzo Salzano for useful suggestions and 
discussion. 


22 


References 


[1] Riess A.G. et al., (Supernova Search Team), Observational Evidence from Supernovae for 
an Accelerating Universe and a Cosmological Constant, Astron. J. 116 (1998) 1009 [arXiv: 
astro-ph/98052011. 

[2] Peiimutter S. et ah, (Supernova Cosmology Project), Measurements of Omega and Lambda 
from 42 High-Redshift Supernovae, Astrophys. J. 517 (1999) 565 [arXiv: astro-ph/9812133 |. 

[3] Eisenstein DJ. et ah. Detection of the Baryon Acoustic Peak in the Large-Scale Corre¬ 
lation Function of SDSS Luminous Red Galaxies, Astrophys. J. 633 (2005) 560 [arXiv: 
astro-ph/05011711; 

Percival WJ. et ah. Measuring the Baryon Acoustic Oscillation scale using the Sloan Digital 
Sky Survey and 2dF Galaxy Redshift Survey, Mon. Not. Roy. Astron. Soc. 381 (2007) 1053 
|arXiv:0705.3323 [astro-ph]]. 

[4] Ade P.A.R. et al. (Planck collaboration), Planck 2013 results. XVI. Cosmological parame¬ 
ters, Astron. Astrophys. 571 (2014) 66 [arXiv: 1303.5076 [astro-ph]]. 

[5] Tegmark M. et ah, (SDSS), Cosmological parameters from SDSS and WMAP, Phys.Rev.D 
69(2004) 103501 [arXiv: astro-ph/0310723]. 

[6] Carroll S.M., The Cosmological Constant, Living Rev. Rel 4 (2001) 1 |astro-ph/0004075]; 
Padmanabhan T., Cosmological constant-the weight of the vacuum, Phys. Rep. 380 (2003) 
235 [arXiv:hep-th/02122901; 

Copeland E.J., Sami M. & Tsujikawa S., Dynamics of Dark Energy, IJMPD 15 (2006) 
17531 arXiv :hep-th/06030571; 

Sahni V. & Starobinsky A., Reconstructing Dark Energy, IJMPD 15 (2006) 2105 
I arXiv:astro-ph/0610026]; 

Erieman J.A., Turner M.S. & Huterer D., Dark Energy and the Accelerating Universe, Ann. 
Rev. Astron. Astrophys. 46 (2008) 385 [arXiv: 0803.0982]; 

Caldwell R.R. & Kamionkowski M., The Physics of Cosmic Acceleration, Ann. Rev. Nucl. 
Part. Sci. 59 (2009) 397 |arXiv:0903.0866]; 

Weinberg D. H. et al.. Observational probes of cosmic acceleration, Phys. Rep. 530 (2013) 
87 |arXiv:1201.2434] 

[7] Riess A. G. et al.. New Hubble Space Telescope Discoveries of Type la Supernovae at z > 1: 
Narrowing Constraints on the Early Behavior of Dark Energy, Astrophys. J. 659 (2007) 98 
[arXiv: astro-ph/0611572| 


23 


[8] Shapiro C. & Turner M.S, What Do We Really Know about Cosmic Acceleration?, Astro- 
phys. J. 649 (2006) 563 |arXiv:astro-ph/0512586]. 

[9] Melchiorri A., Pagano L. & Pandolfi S., When did cosmic acceleration start?, Phys. Rev. D 
76(2007) 041301 |arXiv:0706.1314] 

[10] Rapetti D. et ah, A kinematical approach to dark energy studies, Mon. Not. Roy. Astron. Soc. 
375 (2007) 1510 |arXiv:astro-ph/0605683 | 

[11] Nair R., Jhingan S & Jain D., cosmokinematics: a joint analysis of standard candles, rulers 
and cosmic clocks, JCAP 1201 (2012) 018 |arXiv:1109.4574|. 

[12] Wang F.Y. &: Dai Z.G., Constraining the cosmological parameters and transition redshift 
with gamma-ray bursts and supernovae, Mon. Not. Roy. Astron. Soc. 368 (2006) 371 [arXiv: 
astro-ph/05122791 

[13] Xu Lixin, Li W. & Lu J., Constraints on kinematic model from recent cosmic observations: 
SNla, BAO and observational Hubble data, JCAP 0907 (2009) 031 [arXiv: 0905.4552]. 

[14] Cunha J.V. & Lima J.A.S. Transition redshift: new kinematic constraints from supernovae, 
Mon. Not. Roy. Astron. Soc. 390 (2008) 210 |arXiv:0805.12611. 

[15] Gong Y. & Wang A., Reconstruction of the deceleration parameter and the equation of state 
of dark energy, Phys. Rev. D 75 (2007) 043520 | arXiv: astro-pli/0612196]. 

[16] Lima J.A.S., Holanda R.F.L. & Cunha J.V., Are Galaxy Clusters Suggesting an Accelerating 
Universe Independent ofSNe la and Gravity Metric Theory? (2009) [arXiv: 0905.2628] 

[17] Farooq O., Crandall S. &: Ratra B., Binned Hubble parameter measurements and the cosmo¬ 
logical deceleration-acceleration transition, Phys. Lett. B726 (2013) 72 |arXiv:1305.19571 

[18] Farooq O. &: Ratra B., Hubble Parameter Measurement Constraints on the Cos¬ 
mological Deceleration-Acceleration Transition Redshift, Astrophys. J. 766 (2013) 01 
|arXiv:130L5243| 

[19] Riess A.G. et ah, Type la Supernova Discoveries at z > Ifroni the Hubble Space Telescope: 
Evidence for Past Deceleration and Constraints on Dark Energy Evolution, Astrophys. J. 
607 (2004) 665 |arXiv:astro-ph/0402512] 

[20] Elgaroy O. & Multamaki T, Bayesian analysis of Eriedmannless cosmologies, JCAP 0609 
(2006) 002 [arXiv:astro-ph/0603053] 

[21] Guimaraes A.C.C., Cunha J.V. & Lima J.A.S., Bayesian analysis and constraints on kine¬ 
matic models from union SNIa, JCAP 0910 (2009) 010 | arXiv :0904.35501 


24 


[22] Turner M. & Riess A., Do Type la Supernovae Provide Direct Evidence for Past Decelera¬ 
tion of the Universe?, Astrophys. J. 569 (2002) 18 |arXiv:astro-ph/01060511. 

[23] Cai R.G. &: Tuo Z., Detecting the cosmic acceleration with current data, Phys. Lett. B 706 
(2011) 116 |arXiv:1105.1603] 

[24] Lima J.A.S. et al., Is the transition redshift a new cosmological number? (2012) 
|arXiv:1205.4688| 

[25] Del Campo S. et al., Three thermodynamically based parametrizations of the deceleration 
parameter, Phys. Rev. D 86 (2012) 083509 [arXiv: 1209.3415] 

[26] Vargas dos Santos M., Reis R. R. R &: Waga 1., Constraining cosmic deceleration- 
acceleration transition with type la supernova, BAO/CMB and H(z) data (2015) 
[arXiv:1505.03814] 

[27] Huterer D.& Starkman G., Parametrization of Dark-Energy Properties: A Principal- 
Component Approach, Phys. Rev. Lett. 90 (2003) 031301 [arXiv:astro-ph/0207517 | 
Shafieloo A., Model-independent reconstruction of the expansion history of the Uni¬ 
verse and the properties of dark energy, Mon. Not. Roy. Astron. Soc. 380 (2007) 1573 
I arXiv:astro-ph/0703034] 

Holsclaw T. et al., Nonparametric reconstruction of the dark energy equation of state, Phys. 
Rev. D 82 (2010) 103502 |arXiv: 1009.5443] 

Holsclaw T. et al., Nonparametric Dark Energy Reconstruction from Supernova Data, Phys. 
Rev. Lett. 105 (2010) 241302 [arXivHOl 1.30791 

Lazkoz R., Salzano V. & Sendra 1., Revisiting a model-independent dark energy reconstruc¬ 
tion method, Euro Phys. Jr. C 72 (2012) 2130 [arXiv: 1202.46891 

Seikel M., Clarkson C. & Smith M., Reconstruction of dark energy and expansion dynamics 
using Gaussian processes, JCAP 06 (2012) 036 [arXiv: 1204.2832] 

Shafieloo A., Crossing statistic: reconstructing the expansion history of the universe, JCAP 
08 (2012) 002 [arXiv: 1204:1109] 

Nesseris S. & Garcia-Bellido J., Comparative analysis of model-independent methods for 
exploring the nature of dark energy, Phys. Rev. D 88 (2013) 063521 [arXiv:1306.4885] 
Sapone D., Majerotto E. & Nesseris S., Curvature versus distances: Testing the ELRW 
cosmology, Phys. Rev. D 90 (2014) 023012 [arXiv: 1402.2236] 

Li zhengxiang et ah. Constructing a cosmological model-independent Hubble diagram of 
type la supernovae with cosmic chronometers (2015) [arXiv: 1504.03269 [ 


25 


Vitenti S.D.P. & Penna-Lima M., A general reconstruction of the recent expansion history 
of the universe (2015) |arXiv: 1505.018831 

[28] W. Clewelandjiobust Locally Weighted Regression and Smoothing Scatterplots, J. Amer 
Statist. Assoc. 74 (1979) 82936 

[29] W. Cleveland & S. J. Devlin, Locally Weighted Regression: An Approach to Regression Anal¬ 
ysis by Local Fitting, J. Amer. Statist. Assoc. 83 (1988) 596 

[30] Montiel A. et al, Nonparametric reconstruction of the cosmic expansion with local re¬ 
gression smoothing and simulation extrapolation, Phys. Rev. D 89 (2014) 043007 [arXiv: 
1401.4188] 

[31] Amendola L. & Tsujikawa S., Dark Energy: Theory and Observations (1st edition), Cam¬ 
bridge University Press (2010). 

[32] Huterer D. &: Turner M.S., Prospects for probing the dark energy via supernova distance 
measurements, Phys. Rev. D 60 (1999) 081301 

[33] Efstathiou G., Constraining the equation of state of the Universe from distant Type la Super¬ 
novae and Cosmic Microwave Background Anisotropies, Mon. Not. Roy. Astron. Soc. 342 
(2000)810 

[34] Simon J., Verde L. & Jimenez R., Constraints on the redshift dependence of the dark energy 
potential, Phys. Rev. Dll (2005) 123001 [arXiv:astro-ph/0412269] 

[35] Jimenez R. et al., Synthetic stellar populations: Single stellar populations, stellar inte¬ 
rior models and primordial proto-galaxies, Mon. Not. Roy. Astron. Soc. 349 (2004) 240 
I arXiv:astro-ph/0402271 ] 

[36] McCarthy P.J. et al.. Evolved Galaxies at z > 1.5 from the Gemini Deep Deep Sur¬ 
vey: The Formation Epoch of Massive Stellar Systems, Astrophys. J. Lett. 614 (2004) L9 
[ arXiv:astro-ph/0408367] 

[37] Dunlop J. et al., A 3.5-Gyr-old galaxy at redshift 1.55, Nature 381 (1996) 581 

[38] Nolan L. et al., F stars, metallicity, and the ages of red galaxies at z > 1, Mon. Not. Roy. 
Astron. Soc. 341 (2003) 464 [arXiv:astro-ph/0103450] 

[39] Samushia L. et al.. Constraints on dark energy from the lookback time versus redshift test, 
Phys. Lett. B 693 (2010) 509 [arXiv: 0906.2734] 

[40] Liu G. et al.. The Age-redshift Relation for Luminous Red Galaxies Obtained from Full 
Spectrum Fitting and its Cosmological Implications, Astrophys. J. 758 (2012) 107 [arXiv: 
1208.6502] 


26 


[41] Riess A.G. et al., A 3% Solution: Determination of the Hubble Constant with the Hubble 
Space Telescope and Wide Field Camera 3, Astrophys. J. 730 (2011) 119 |arXiv: 1103.2976 j 

[42] Lima J.A.S. & Cunha J.V., A 3% Determination Of Hq at intermediate redshifts, Astrophys. 
J. 781 (2014) L38 [arXiv: 1206.0332] 

[43] Ade R A.R. et al. (Planck collaboration), Planck 2015 results. XIII. Cosmological parameters 
[arXiv:1502.01589] 

[44] Narayan R. &: Bartelmann M., Lectures on Gravitational Leasing (1996) 
[ arXiv: astro-ph/9606001 ] 

[45] Cao S. et al., Constraints on cosmological models from strong gravitational leasing systems, 
JCAP 1203 (2012) 016 [arXiv: 1105.6226] 

[46] Ofek E.O. et al.. The redshift distribution of gravitational lenses revisited: con¬ 
straints on galaxy mass evolution, Mon. Not. Roy. Astron. Soc. 343 (2003) 639 
]arXiv:astro-ph/0305201 ] 

[47] Liao K & Zhu Z-H, Constraints onf(R) cosmologies from strong gravitational leasing sys¬ 
tems, Phys. Lett. B 714 (2012) 1 ]arXiv:1207.2552 ] 

[48] Chen Y. et al.. Constraints on a phiCDM model from strong gravitational leasing and up¬ 
dated Hubble parameter measurements, JCAP 1502 (2015) 010 [arXiv: 1312.1443] 

[49] Betoule M. et al.. Improved cosmological constraints from a joint analysis of the SDSS-II 
and SNLS supernova samples, Astron. Astrophys. 568 (2014) 32 ]arXiv:140L4064 ] 

[50] Li Zhengxiang, Ding X. & Zhu Zong-Hong Unbiased constraints on the dumpiness of the 
Universe from standard candles, Phys. Rev. D 91 (2015) 083010 [arXiv: 1504.03482] 

[51] Liddle A. R., How many cosmological parameters?, Mon. Not. Roy. Astron. Soc. 351 (2004) 
L49 ]arXiv:astro-ph/0401198 ] 

[52] Davis T.M. et al.. Scrutinizing Exotic Cosmological Models Using ESSENCE Super¬ 
nova Data Combined with Other Cosmological Probes, Astrophys. J. 666 (2007) 716 
] arXiv:astro-ph/0701510] 

[53] Szydlowski M. et al., AIC, BIC, Bayesian evidence against the interacting dark energy 
model, Euro Phys. Jr. C 75 (2015) 1 ]arXiv:080L0638 ] 

[54] Dantas M.A. et al.. Time and distance constraints on accelerating cosmological models, 
Phys. Lett. B 699 (2011) 239 [arXiv: 1010.0995 ] 


27 


[55] Wei Jun Jie et al., The Age—Redshift Relationship of Old Passive Galaxies, Astron. J. 150 
(2015) 13 |arXiv:1505.07671| 

[56] Biesiada M., A. Piorkowska & B. Malec, Cosmic equation of state from strong gravitational 
leasing systems, Mon. Not. Roy. Astron. Soc. 406 (2010) 1055 |arXiv:l 105.09461 

[57] Melia R, Wei Jun-Jie & Wu Xue-Feng, A Comparison of Cosmological Models Using Strong 
Gravitational Tensing Galaxies, Astron. J. 149 (2015) 1 [arXiv: 1410.0875 | 

[58] Yuan C.C. & Wang F.Y., Cosmological test using strong gravitational leasing systems, Mon. 
Not. Roy. Astron. Soc. 452 (2015) 2423 [arXiv: 1506.08486] 

[59] Treu T. et al., The Sloan Lens ACS Survey. II. Stellar Populations and Internal Structure of 
Early-Type Lens Galaxies, Astrophys. J. 640 (2006) 662 [arXiv:astro-ph/0512044] 

[60] Wang N. Sz Xu L., Strong Gravitational Tensing and its Cosmic Constraints, Mod. Phy. Lett. 
A 28 (2013) 1350057 [arXiv: 1309.0887] 

[61] Van de Ven G., Van Dokkiim P.G. & Franx M. The Fundamental Plane and the evolution of 
the M/L ratio of early-type field galaxies up to z^l Mon. Not. Roy. Astron. Soc. 344 (2003) 
924 [arXiv: astro-ph/0211566] 

[62] Biesiada M. et al., Dark energy constraints from joint analysis of standard rulers and stan¬ 
dard candles. Res. Astron. Astrophys. 11 (2011) 641 

[63] Dalai N., Hennawi J.F & Bode P, Noise in Strong Tensing Cosmography, Astrophys. J. 621 
(2005) 99 [arXiv: astro-ph/0409028] 

[64] Bar-kana R., Effect of Large-Scale Structure on Multiply Imaged Sources, Astrophys. J. 468 
(1996) 17 [arXiv: astro-ph/9511056] 

[65] Keeton C.R., Kochanek C.S. & Seljak U., Shear and Ellipticity in Gravitational Lenses, 
Astrophys. J. 482 (1997) 604 [arXiv: astro-ph/9610163] 

[66] Grillo C., Lombardi M. & Berlin G., Cosmological parameters from strong gravitational 
leasing and stellar dynamics in elliptical galaxies, Astron. Astrophys. 477 (2008) 397 
[arXiv:0711.08821 

[67] Sandage A.R., Cosmology:A search for two numbers, Phys. Today 23 (1970) 34 

[68] Neben A.R. & Turner M.S., Beyond Hq and qq: Cosmology is no longer just two numbers, 
Astrophys. J. 769 (2013) 133 

[69] Cardenas V. H., Bernal C. & Bonilla A., Cosmic slowing down of acceleration using fgas, 
Mon. Not. Roy. Astron. Soc. 433 (2013) 3534 [arXiv:1306.0779 j 


28 


[70] Margana J., Cardenas V.H. & Motta V., Cosmic slowing down of acceleration for several 
dark energy paranietrizations, JCAP 1410 (2014) 017 |arXiv: 1407.1632] 

[71] Schwab J., Boltan A.S. & Rappaport S.A., Galaxy-Scale Strong-Lensing Tests of Gravity and 
Geometric Cosmology: Constraints and Systematic Limitations, Astrophys. J. 708 (2010) 
750 [arXiv:0907.49921 

[72] Rasanen S., Bolejko K. & Finoguenov A., A new test of the FLRW metric using distance 
sum rule, Phys. Rev. Lett. 115 (2015) 101301 [arXiv: 1412.4976 j 

[73] Serjeant S., Up to 100,000 reliable strong gravitational lenses in future dark energy experi¬ 
ments, Astrophys. J. 793 (2014) LIO 

[74] www.sdss3.org/surveys/boss.php 


29 


